function R1 = findR2(Q1,i_new,R,n)
    global I fL1
    %FL1 = @function_F;
    %fL1 = @function_Poiss;
    %载入R矩阵
    N=size(R,2);
    R1=[R(1:max(i_new-1,1),:);zeros(I-max(i_new-1,1),N)];
    
    %如果新的病人安排到slot1，重新计算R(1,:)，并将i_new修改为2
    if i_new==1
        for k = 1:n 
            R1(1,k+1)=sum(fL1(1:n-k+1).*Q1(1,k+1:n+1));
        end
        R1(1,1)=1-sum(R1(1,2:n+1));%利用i时期R概率和为1计算R1(1,1);
        i_new=i_new+1;%将i_new的值更改为2。
    end
    
    %更新R的i_new到I行。
    for i=i_new:I  
    %计算R(i,:)
    for k=1:n
        for L=0:n-k%此处L指L的值
            for m=0:k+L
            R1(i,k+1)=R1(i,k+1)+fL1(L+1)*Q1(i,m+1)*R1(i-1,k+L-m+1);
            end
        end
    end
    R1(i,1)=1-sum(R1(i,2:n+1));%利用i时期R概率和为1计算R1(i,1);
    end
return;